Universal spectral statistics of Andreev billiards: semiclassical approach 
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The symmetry classification of complex quantum systems has recently been extended beyond 
the Wigner-Dyson classes. Several of the novel symmetry classes can be discussed naturally in the 
context of superconducting-normal hybrid systems such as Andreev billiards and graphs. In this 
paper, we give a semiclassical interpretation of their universal spectral form factors in the ergodic 
limit. 
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I. INTRODUCTION 

Based on early work of Wigner 1], Dyson proposed 
a classification of complex quantum systems according 
to their behavior under time reversal and spin rotations. 
The ergodic limits of the proposed symmetry classes are 
described by the Gaussian orthogonal, unitary, and sym- 
plectic ensembles (GOE, GUE, and GSE) of random- 
matrix theory. These were initially motivated by atomic 
nuclei and have since been applied successfully to a large 
variety of systems^most notably chaotic and disordered 
quantum systems 'j] . More recently, an additional seven 
symmetry classes have been identified ^ , which are nat- 
urally realized in part by Dirac fermions in random gauge 
fields (chiral classes) Q and in part by quasiparticles in 
disordered mesoscopic superconductors or supercon- 
ducting- normalconducting (SN) hybrid systems The 
common new feature of the new symmetry classes is a 
mirror symmetry in the spectrum: if E is in the spec- 
trum, so is ~E. The corresponding Gaussian random- 
matrix ensembles differ from the Wigner-Dyson ensem- 
bles in so far as their spectral statistics, while still uni- 
versal, is no longer stationary under shifts of the energy 
due to additional discrete symmetries. 

Much insight into the range of validity of the Wigner- 
Dyson random-matrix ensembles has been gained from 
the semiclassical approach to the spectral statistics of 
chaotic quantum systems, based on Gutzwiller's trace 
formula In a seminal paper Ifl], Berry gave a semi- 
classical derivation of the spectral form factor of chaotic 
quantum systems for the Wigner-Dyson ensembles, par- 
tially reproducing the results of random-matrix theory 
and clarifying its limitations. In this paper we provide 
such a semiclassical interpretation, based on Andreev 
systems, of a generalized form factor for the Gaussian 
random-matrix ensembles associated with the new sym- 
metry classes termed C and CI (the pertinent Gaussian 
ensembles will be referred to as C-GE and CTGE). 
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There have been several attempts to a pply semiclassi- 
cal theory to SN hybrid systems 0, [3| ■ Melsen et 
al. [1^ pointed out that the gap induced by the proxim- 
ity effect in a billiard coupled to a superconducting lead 
(Andreev billiard - cf. Fig. lb) is sensitive to whether the 
classical dynamics of the (normal) billiard is integrable 
or chaotic. These authors showed that the proximity- 
induced hard gap in the chaotic case is not fully re- 
produced by semiclassical theory, the reasons for which 
have been discussed further in Ref. j^J. All these sys- 
tems have in common that due to the presence of the 
superconductor the combined electron-hole dynamics is 
no longer chaotic even if the corresponding normal (non- 
superconducting) billiard is chaotic By contrast, 
we focus here on a semiclassical approach to SN hybrid 
systems where the combined electron-hole dynamics re- 
mains chaotic even in the presence of the superconductor. 
Such systems exhibit the universal spectral statistics of 
the Gaussian random-matrix ensembles for the new sym- 
metry classes. We identify the class of periodic orbits 
contributing to universal features in the density of states 
near the Fermi energy and show that in this case, semi- 
classics reproduces the spectral statistics predicted by 
random-matrix theory. 



II. UNIVERSAL SPECTRAL STATISTICS 

We briefly summarize the pertinent random-matrix re- 
sults for the Gaussian ensembles corresponding to the 
new symmetry classes. For the Wigner-Dyson ensem- 
bles (GUE, GOE, GSE), the average density of states 
is nonuniversal and random-matrix theory makes univer- 
sal predictions only about spectral fluctuations in the 
ergodic limit such as the correlation function 
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of the deviations Sp{E) of the density of states p{E) from 
its mean value {p{E)). A central quantity is the spectral 
form factor 
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The ergodic limit of the new symmetry classes differs 
from the Wigner-Dyson case by the fact that even the 
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average density of states has universal features close to 
the Fermi energy /x. Thus, in this case, we define a gen- 
eralized spectral form factor by the Fourier transform of 
the expectation value of the (oscillating part of the) den- 
sity of states: 
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where E is the energy (measured relative to /i). For the 
ensemble C-GE (class C is invariant under spin rotations, 
while time reversal is broken) , this form factor is 
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Here, tn = 27r?ipav is the Heisenberg time defined in 
terms of the mean density of states pav sufliciently far 
from the Fermi energy (the oscillating part of the den- 
sity of states is defined as Sp{E) = p{E) — ps_y). Semiclas- 
sically pav corresponds to Weyl's law. For the ensemble 
CI-GE (class CI differs from C by invariance under time- 
reversal), the short-time expansion is |2|| 
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Wigner-Dyson statistics can be applied even to a sin- 
gle chaotic system by exploiting a spectral average. By 
contrast, the new symmetry classes require an ensemble 
average since they have universal features in the vicinity 
of special energies (Fermi energy p). For billiards one 
may average over shapes. 

Before entering into the semiclassical analysis for the 
new symmetry classes, we briefly review the semiclassical 
derivation of the usual spectral form factor of the GUE. 
There one starts from the Gutzwiller trace formula, that 
relates the oscillatory contribution 6p{E) to the density 
of states to a sum over periodic orbits p, 



Sp{E) = ^Rc V tpApc'^" 
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Here, 



Sp denotes the classical action of the orbit, Ap de- 



notes its stability amplitude, and tp is the primitive orbit 
traversal time. The explicit factor tp arises because the 
traversal of the periodic orbit can start anywhere along 
the orbit. Inserting this expression into the definition 
of the spectral form factor, and employing the diagonal 
approximation, one finds 
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Finally averaging over some time interval At and 
using the Hannay-Ozorio-de-Almeida sum rule 
^tpe[t,t+At] |^^|2 ^ ^^^^ Qj^g obtains the resuh 
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FIG. 1; (a) Andreev scattering with a perpendicular mag- 
netic field, (b) Andreev billiard - a part of the boundary 
is connected to a superconductor where Andreev reflection 
(retroflection) takes place. An example of a self-dual orbit is 
shown - the SN interface is hit three times by an incoming 
electron and three times by an incoming hole. Each part of 
the trajectory is traversed twice in the same direction - once 
as an electron and once as a hole. The singly traversed orbit 
without electron-hole labels is a periodic orbit of the virtual 
billiard and hits the SN interface three times. 



vahd for to <^ t tn, where to is the period of the 
shortest periodic orbit. This result agrees with the short- 
time behavior of the spectral form factor predicted by the 
GUE. 



III. SEMICLASSICAL APPROACH TO 
MAGNETIC ANDREEV BILLIARDS 

We now turn to Andreev billiards - the central theme 
of this paper. The novel element in SN hybrid systems 
is Andreev reflection converting electrons into holes (and 
vice versa) at the interface to the superconductor (see 
Fig. la). In this process, the incoming electron (hole) 
acquires a phase — ie~'" (— ie'"), where a is the phase 
of the superconducting order parameter A T^. In the 
absence of a magnetic field, electrons (holes) sufficiently 
close to the Fermi energy (i? <C |A| <^ p,) are reflected as 
holes (electrons) which then retrace the electron (hole) 
trajectory backwards (retroflection). In chaotic billiards, 
essentially all trajectories eventually hit any given part of 
the boundary. Thus, if the billiard is then coupled to a su- 
perconductor any quasiparticle eventually hits the super- 
conducting interface, leading to a periodic orbit bouncing 
back and forth between two points on the superconduct- 
ing interface. It follows that a conventional chaotic bil- 
liard (without magnetic field) that is coupled to a super- 
conductor has a combined electron-hole dynamics that 
is no longer chaotic. Instead, the resulting trajectories 
are all periodic, leading to nonuniversal behavior such as 
the proximity- induced hard gap 0, for time-reversal 
invariant systems. 

One expects to recover universal spectral statistics only 
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if the combined electron-hole dynamics is chaotic and pe- 
riodic orbits are isolated as in conventional chaotic (hy- 
perbolic) systems. In Andreev billiards this occurs nat- 
urally when time-reversal symmetry is broken by a per- 
pendicular magnetic field (symmetry class C). In this 
case the retroflected hole (electron) does not retrace the 
trajectory of the incoming electron (hole), as electron 
and hole trajectories are curved in the same direction 
(cf. Fig. la). This allows one to express the density of 
states semiclassically by a Gutzwiller-type trace formula 
as a sum over the isolated periodic orbits of the combined 
electron-hole dynamics of the Andreev billiard, 

Sp{E) = ^ReV tpApe'^''(^)/'"'+'x. (9) 
nn, ^ — ' 

p 

The orbit amplitudes Ap are products of electron and 
hole contributions, 

Ap = A(;^A'^\ (10) 

while the orbit actions are sums of electron and hole ac- 
tions, 

Sp{E)^Si^HE) + Sl!^\E). (11) 

The factor tp again reflects the arbitrary starting point of 
the orbit and x denotes the accumulated Andreev phases. 

Coherent contributions to the form factor can be ex- 
pected from the periodic orbits that retrace the same tra- 
jectory in the same direction with the roles of electrons 
and holes interchanged. Such self-dual orbits are invari- 
ant under electron-hole conjugation and the dynamical 
contributions to their action largely cancel due to the 
relation sI''\e) = -S^''\-E), so that Sp{E) ~ Etp. 
Moreover, the amplitudes of electron and hole are just 
complex conjugates of one another, giving Ap = \A^p'^\'^. 
The accumulated Andreev phase is (— i)^* = —1 with s 
an odd integer. Keeping only the self-dual periodic orbits 
- the self-dual approximation ~ we find 

Sp{EU = -4rRe Vtpl^'^f e'^*-/^ (12) 
p 

For the generalized form factor, this leads to 

K{tU = -2j2tp\Al;''>\^S(t~tp). (13) 

p:sd 

This expression reveals the similarity to the diagonal ap- 
proximation for the Wigner-Dyson form factors. How- 
ever, only one factor tp arises. 

The Hannay-Ozorio-de- Almeida sum rule does not ap- 
ply directly to the sum over self-dual orbits (being a 
sum over amplitudes of a subclass of periodic orbits). 
To deal with this difficulty, we introduce a virtual bil- 
liard with the same dynamics as the Andreev billiard 
except that there is no particle-hole conversion at the 
SN interface. Thus, the virtual billiard is an ordinary 



chaotic billiard with unusual reflection conditions at the 
SN interface (retroflections). Primitive periodic orbits of 
the virtual billiard involve either even or odd numbers 
of retroflections. Reintroducing electron-hole conversion, 
one observes that even orbits lead to non-self-dual peri- 
odic orbits in the Andreev billiard. By contrast, twofold 
traversals of odd orbits are periodic and self dual in the 
Andreev billiard as the roles of electron and hole are in- 
terchanged in the second traversal (see Fig. lb). We can 
now interpret the sum over self-dual orbits in Eq. l(T^ 
as a sum over odd orbits of the virtual billiard. Since on 
average half of its orbits are odd, the Hannay-Ozorio-de- 
Almeida sum rule for the virtual billiard gives 

tpe[t,t+At] 

p:sd 

Thus 

K{tU = -1 (15) 

in agreement with the random-matrix result predicted by 
C-GE (0J for short times. The self-dual approximation 
is expected to hold for to,tA <^ t <^ tn where to is the 
traversal time of the shortest periodic orbit and the 
Andreev time (typical time until electron-hole conversion 
takes place). 

IV. SPECTRAL STATISTICS FOR ANDREEV 
GRAPHS 

The semiclassical calculation of the form factor be- 
comes particularly transparent and explicit for quantum 
graphs which were recently introduced |l6j| as simple 
quantum chaotic systems. Introducing Andreev reflec- 
tion as a new ingredient, we show semiclassically that the 
form factor of the resulting Andreev graph takes on the 
universal result. A quantum graph consists of vertices 
connected by bonds. A particle (electron/hole) propa- 
gates freely on a bond and is scattered at a vertex accord- 
ing to a prescribed scattering matrix. For definiteness, 
we discuss star graphs with TV bonds of equal length L. 
These have one central vertex and N peripheral vertices. 
Each bond connects the central vertex to one peripheral 
vertex (cf. Fig. 2). 

Andreev (star) graphs are obtained by introducing 
(complete) electron-hole conversions at the peripheral 
vertices, while the central vertex preserves the particle 
type. The quantization condition is 

det(5(fc) - 1) = 0, (16) 

with the unitary N x N matrix 

S{k) = ScCD^CS*cCD+C. (17) 

Here Sc (Sq) is the central scattering matrix for an elec- 
tron (hole). For definiteness, we choose ^3 

ScM = -^e'-'^'/^, (18) 
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FIG. 2: Star graph with five peripheral vertices connected to 
superconductors. 



where Sc by itself does not break time-reversal symme- 
try. The matrix 



C 



(19) 



contains the phases accumulated when the quasiparticle 
propagates along the bonds {k is the wave number mea- 
sured from the Fermi wave number). Finally, 



D 



± 



-idiag(e=Fi"-) 



(20) 



contains the Andreev phases accumulated at the vertices, 
where denotes the order-parameter phase at peripheral 
vertex i. Time-reversal symmetry is obeyed if all Andreev 
phases are either ai — or a.i ~ tt, but is broken oth- 
erwise. Accordingly, we build ensembles corresponding 
to the symmetry classes C (uncorrelated Andreev phases 
ai with uniform distributions in the interval [0,27r)) and 
CI (uncorrelated Andreev phases taking values ai = 
or ai = TT with equal probability). Numerically com- 
puted ensemble averages are in excellent agreement with 
random-matrix results from C-GE and CI-GE, as shown 
in Fig. 3. 

Following previous work on quantum graphs |l6j |. we 
write the density of states in k space as 



p{k) = pav + Sp{k) 



(21) 



with pav = 2NL/tt and obtain the exact trace formula 
Sp{k) = -Ke y tpApe'^^+'^ (22) 



as a semiclassical sum over periodic orbits p of the 
graph. Here, periodic orbits are defined as a sequence 
11,12, ■■■ of peripheral vertices, with cyclic permuta- 
tions identified. Since the particle type changes at the 
peripheral vertices, the sequences must have even length 
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FIG. 3: Form factors for class C (top) and CI (bottom) cal- 
culated numerically for a star graph (full lines) with A'" = 100 
bonds (averaged over 50 000 realizations and a short time in- 
terval At '^tn) and as obtained from the Gaussian random- 
matrix ensembles C-GE and CI-GE (dashed lines) in dimen- 
sionless time r = t/tH- The insets give the coefficients Km as 
a function of "time" m/N . 



I ~ 2m. The primitive traversal "time" |0 of a periodic 
orbit is tp — AmL/r (where r is the repetition number), 
the stability amplitude is Ap = 1/7V™, and the action is 



2m 



+lo_^j'j+l 



Sp = AmkL + ^(-l)^+i27r 
The accumulated Andreev phase is 



X = -rmr - ^(-1)^+^0^^ . 
Then, the form factor becomes 



N 



(23) 



(24) 



K{t) = 2 / 

J — C 



dke~'^' 



{5p{k)) 



m— 1 



(25) 



with the Heisenberg time tn = iLN ((•) denotes the 
average over Andreev phases). The coefficients can be 
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FIG. 4: Periodic orbits contributing in the self-dual approxi- 
mation (at m — 5). The vertices in the diagram correspond to 
peripheral vertices of the original star graph, full and dashed 
lines represent electron and hole propagation. In class C, only 
the left diagram contributes. In class CI, the right diagram 
gives m additional contributions as the turning point can be 
any of the m vertices. 



written as a sum over periodic orbits Pm with 2m An- 
dreev reflections: 

X,„-2 5]^(V'''^'^°^+''')- (26) 

Km can be viewed as a form factor in discrete time m/N. 

For graphs in class C, only those periodic orbits sur- 
vive the average over Andreev phases that visit each pe- 
ripheral vertex an even number of times - half as incom- 
ing electron and half as incoming hole. In the self-dual 
approximation, only those orbits contribute whose total 
phase due to the scattering matrix of the central vertex 
vanishes. As the phase factors due to scattering between 
bonds i and j for electrons and holes are complex con- 
jugates of one another, this requires that the periodic 
orbits contain equal numbers of scatterings from i to j 
as electron and hole. This leads to the orbits sketched in 
Fig. 4: An odd number of peripheral vertices are visited 
twice, once as an electron and once as a hole. First, the 
peripheral vertices are visited once, alternating between 
electrons and holes, and subsequently the vertices are 
visited again in the same order but with the roles of elec- 
trons and holes interchanged. Thus, these orbits have the 
same structure as the self-dual orbits discussed above for 
the Andreev biUiard. We have Ap = 1/iV™, Sp = AmkL, 
and X = mTT. The number of such orbits of length 2m is 
N"^ /m, where the denominator m reflects the identifica- 
tion of cyclic permutations of peripheral vertices. With 
these ingredients, we find the short-time result 

i^,'^,.d--l + (-iy" ^fd(i)--l, (27) 

where K(t) is the time-averaged form factor. This repro- 
duces the result predicted by C-GE. 

For class CI, the average over Andreev phases requires 
only an even number of visits to each vertex. In the self- 
dual approximation, this leads to additional orbits (see 
Fig. 4) and to the result 

KgU = -1 + (-l)'"(2m + 1) ^ K^ht) = -1, (28) 

the leading order term for short times predicted by the 
corresponding random-matrix ensemble CTGE. 



V. ANDREEV BILLIARDS WITHOUT 
MAGNETIC FIELD 

Finally, we come back to Andreev billiards without 
magnetic field (class CI). As explained above (in section 
imp holes necessarily retrace the electron trajectory, thus 
leading to non-isolated periodic orbits and nonuniversal 
spectral statistics (hard gap). Universal spectral statis- 
tics of CI can, however, be found in such Andreev bil- 
liards with N one-channel leads. The reason for this is 
that Andreev billiards with N leads containing one chan- 
nel each can be mapped to star graphs. The quantization 
condition for Andreev billiards with N leads is ^3 

det{S{E) - 1) = 0. (29) 

Here S{E) is the N x N Andreev billiard scattering ma- 
trix 

S{E) = SNc{E)D^S%ci-E)D+, (30) 

with SNciE) the scattering matrix describing the cou- 
pling of the N channels by the normal region. The matri- 
ces D± describing the Andreev scattering in the leads are 
diagonal, D± = — idiag(e^'"'), with a specific Andreev 
phase ai for each lead. Time reversal invariance demands 
= or = tt. Then a detailed correspondence be- 
tween billiard and star graph is obtained by substituting 
CScC Snc{E) and CS^C ^ S'iVc(-^) (with a more 
general central scattering matrix). Thus, the form factor 
of these billiards can be obtained in the self-dual approx- 
imation in complete analogy with the star graph. 

VI. CONCLUSIONS 

We considered the universal spectral statistics for er- 
godic SN hybrid systems belonging to the new symmetry 
classes, in the semiclassical approximation. While it was 
known that semiclassics has problems in some types of 
Andreev systems ^3 we showed both for billiards 
and for quantum graphs that the universal spectral statis- 
tics of the random-matrix ensembles C-GE and CTGE as 
reflected by the appropriately generalized form factor is 
correctly reproduced by semiclassical theory. An impor- 
tant condition for finding the universal spectral statistics 
is that the combined electron-hole dynamics of the An- 
dreev system is classically chaotic. In particular, this re- 
quires that the hole does not retrace the trajectory of the 
incoming electron. In class C, this is naturally the case in 
magnetic Andreev billiards. We related the universal fea- 
tures in the density of states to self-dual periodic orbits 
which are invariant under electron-hole exchange. Our 
results clarify under which conditions to expect spectral 
statistics described by the novel random-matrix ensem- 
bles. 

The results presented can be extended to the symmetry 
classes D,Dlll and the chiral classes. We also note that 
our results for Andreev graphs remain valid for a rather 
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large class of central scattering matrices Sc- Finally, 
by going beyond the self-dual approximation in Andreev 
graphs, it is possible to extract the orbits contributing to 
the form factor to linear order in t (weak localization cor- 
rections). These extensions will be discussed elsewhere 



Acknowledgments 

We thank R. Klesse, S. Nonnenmacher, H. Schanz, U. 
Smilansky, and especially P. Brouwer for helpful discus- 
sions. This work was supported in part by SFB 290. 



[1] E.P. Wigner, Ann. Math. 67, 325 (1958). 

[2] F.J. Dyson, J. Math. Phys. 3, 140, 1199 (1962). 

[3] For a review, see T. Guhr, A. Muller-Groeling, and 

H.A. Weidenmuller, Phys. Rep. 299, 189 (1998). 
[4] M.R. Zirnbauer, J. Math. Phys. 37, 4989 (1996). 
[5] J.J.M. Verbaarschot and L Zahed, Phys. Rev. Lett. 70, 

3852 (1993); J.J.M. Verbaarschot, tbid 72, 2531 (1994). 
[6] A. Altland, B.D. Simons, and M.R. Zirnbauer, 

Phys. Rep. 359, 283 (2002). 
[7] A. Ahland and M.R. Zirnbauer, Phys. Rev. B 55, 1142 

(1997). 

[8] M.C. Gutzwiller, J. Math. Phys. 12, 343 (1971) 
[9] M.V. Berry, Proc. R. Soc. London, Ser. A 400, 229 
(1985). 

[10] J. Melsen, P. Brouwer, K. Frahm, and C. Beenakker, 

Europhys. Lett. 35, 7 (1996). 
[11] D. Taras-Semchuk and A. Altland, Phys. Rev. B 64, 

014512 (2001). 

[12] W. Ihra, M. Leadbeater, J.L. Vega, and K. Richter, 



Eur. Phys. J. B 21, 425 (2001). 
[13] I. Kosztin, D.L. Maslov, and P.M. Goldbart, 

Phys. Rev. Lett 75, 1735 (1995). 
[14] J.H. Hannay and A.M. Ozorio de Almeida, J. Phys. A 

17, 3429 (1984). 
[15] Strictly speaking, the phase of the order parameter is 

not a gauge invariant quantity. For class CI, we assume 

B = with the gauge A = throughout this paper. 

The discussion of class C is valid in any fixed gauge. 

All physical properties such as the form factor are easily 

shown to be gauge independent. 
[16] T. Kottos and U. Smilansky, Phys. Rev. Lett. 79, 4794 

(1997); Ann. Phys. 274, 76 (1999). 
[17] G. Tanner, J. Phys. A 34, 8485 (2001). 
[18] Although the bond length is a dimensionless quantity, we 

nevertheless refer to it loosely as time. 
[19] S. Gnutzmann, B. Seif, F. von Oppen, and M.R. Zirn- 
bauer (unpubhshed). 



